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Abstract. Given a quadratic map Q : K n -» K k defined over a 
computable subring D of a real closed field K, and p G D\Y\, . . . ,Y^] 
of degree d we consider the zero set Z = Z(p(Q(X)),K n ) C K n of 
. . . , X n )) G Z?[-Xi, . . . , X n ]. We present a procedure that com- 
putes, in (dn)°( k ^ arithmetic operations in D, a set S of (real univariate 
representations of) sampling points in K n that intersects nontrivially 
each connected component of Z. As soon as k = o(n), this is faster 
than the standard methods that all have exponential dependence on n 
in the complexity. In particular, our procedure is polynomial-time for 
constant k. In contrast, the best previously known procedure is only 
capable of deciding in n olyk > operations the nonemptiness (rather than 
constructing sampling points) of the set Z in the case of p(Y) = ^ Y- 2 
and homogeneous Q. 

A by-product of our procedure is a bound (dn)°( k ' on the number of 
connected components of Z. 

The procedure consists of exact symbolic computations in D and outputs 
vectors of algebraic numbers. It involves extending K by infinitesimals 
and subsequent limit computation by a novel procedure that utilizes 
knowledge of an explicit isomorphism between real algebraic sets. 
Keywords. symbolic computation, complexity, semialgebraic set, 
quadratic map, univariate representation, infinitesimal deformation. 
Subject classification. 68W30, 13P10, 14Q20, 14Pxx 



1. Introduction and the results 

The algorithmic problem of finding points in real algebraic sets has received 
considerable attention, in particular as it forms a building block for a lot of 
procedures in real algebraic geometry |5j. Even if the algebraic sets one is in- 
terested in are subsets of W 1 , the algorithms with the best known complexity 
bounds use transcendental infinitesimals extending R to perform necessary ge- 
ometric deformations of the sets. Hence it is natural to describe the procedures 
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as operating over an arbitrary real closed field K, with the input data, i.e. the 
polynomials, lying in D[Xi, . . . ,X n ] = D[X], with D C K a computable (in 
sense discussed e.g. in j5J Sect. 8.1]) subring of K. In the case K = M one 
usually assumes D — Z. 

Let 5j's be the connected components of the real algebraic set S = Z(f, K n ), 
/ G D[Jf]. In general the number of the <Sj's is bounded by d 0(jl \ where 
d = degf, and this bound is sharp, see Theorem 7.23, Remark 7.22] and 
Remark 11.71 below. We are interested in point-finding (also called sampling) 
algorithms that produce a finite set of points that intersects each <Sj. Such al- 
gorithms with the best known complexity bounds need at most d°^ arithmetic 
operations in D. Often such sets are exactly what is needed in applications. 
We describe here a procedure that finds a point in each connected component 
of our class of algebraic sets, namely the sets of the form Z(p(Q(X)), K n ), for 
p G D\Yi, . . . , Y k ] of degree degp = d and Q = (Qi(X), . . . , Q k (X)) G D[X] k a 
quadratic map, i.e. degQj < 2 for 1 < j < k with the complexity (dn)°^ h \ 

The result of [3J, that bounded, in particular, the sum of the Betti num- 
bers of the set of real solutions of a system of quadratic equations Q\(X) = 
■■■ = Qk{X) = (that can obviously be written as Z(p(Q(X)), W" 1 ) with 
p(Y) = ^jYf) by a polynomial in k and n of degree 0(k) was perhaps the 
earliest indication that in the case Z(p(Q(X)), K n ) the number of connected 
components has only polynomial dependence on n. However, until the present 
work, an algorithmic procedure with the similar complexity bound (dn)°( k > 
for finding points in Z(p(Q(X)), K n ) was unknown. Even procedures that de- 
cide non-emptiness of Z(p(Q(X)), K n ) in time polynomial in n and d for fixed 
k were, for general p, unknown; in [2] such a procedure was described for 
p(Y) = J2j Yj ', homogeneous Qi's and K = R. 

The technique we use is that of symbolic computation. All the data is 
represented exactly, as (real) algebraic numbers, if necessary. More precisely, 
elements of K n that we compute with are given by real univariate representa- 
tions. The latter are defined as follows. A sign condition for a set of poly- 
nomials V = {Pi,...,P s } C K[Y] is specified by a G {—1,0, 1} S so that 
a = (sign Pi (y), . . . , signP s (F)). Thorn encoding jSJ Lemma 2.38, Sect. 10.4] 
of a root a G IK of / G K[T] is a sign condition a Q on the derivatives of /, 
that is a a = (sign f'(a), . . . , sign /( deg ^ -1 )(a)). Note that a a and / determine 
a G K. Let K denote the algebraic closure of K. A univariate representation of 
u G K is an (m + 2)-tuple 

(1-1) u(T) = (f,go,gi,...,g m ) 

of univariate polynomials in D[T] satisfying u = -4-r(gri(a), . . . , g m (a)) for a 
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root a G K of /, and such that / and go are coprime. Obviously, each u(T) 
can represent as many as deg(/) distinct elements of K . A real univariate 
representation of u G K m is a pair u(T), a a , where a a is the Thorn encoding of 
a root a G K of /. 

The main result of the paper is as follows. 

Theorem 1.2. Let Q = (Qi, • • • , Qk) G D[X X ,..., X n ] k be a quadratic map 
X \— > Q(X), and let p G D\Yi, . . . , Y^} satisfy degp < d. A set of real univariate 
representations u(T), o a of a set of points in Z = Z(p(Q(X)),K n ) meeting each 
connected component of Z can be computed in (dn)°^ arithmetic operations 
in D. The degrees of polynomials in u(T) are bounded by (dn)°^ k \ When 
D — Z, the coefficients of u(T) and the intermediate polynomial data will be 
bounded by (dn)°^ times the bitsize of the input data p, Q. 

From now on whenever we talk about finding points in K n , they are meant 
to be given as real univariate representations. 

Remark 1.3. When IK is archimedean, e.g. K = K, the approximations, 
in the ring of fractions of D, of the point in K n given by a real univariate 
representation can be found efficiently as long as approximations a of a can be 
computed efficiently (indeed, then one can just compute u{6t) ). For instance 
when D = Z one can find an interval X = [a_, a + ] 3 a with a± G Q so that 
a is the only root of f in X, see e.g. JSJ Sect. 10.2]. Once X is known, one 
can compute its repeated (rational) bisections to obtain approximations of a 
of needed precision; the complexity of the latter is analyzed e.g. in 'All (see 
also HQ). 

Note that by connected (component of) semialgebraic set, we mean semial- 
gebraically connected, that is, connected in the semialgebraic topology, (com- 
ponent of) semialgebraic set, see e.g. jHj- It is well-known that for the semial- 
gebraic sets over R semialgebraic connectedness implies connectedness (in the 
usual Euclidean topology), see e.g. Thm. 5.21]. 

Theorem II .21 is proved in Section El by exhibiting a procedure that does the 
claimed task. It immediately implies the following. 

Corollary 1.4. The number of connected components of the set Z is at most 
(dn)°^. □ 

An extra argument, to be published elsewhere, allowed us to show that the 
latter bound holds for the sum of Betti numbers of Z, and not only for the 0-th 
one, i.e. the number of components. As well, one can modify the procedure of 
Theorem 11.21 to prove 
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Theorem 1.5. Under the assumptions of Theorem U.'A computing the exact 
minimum and a minimizer (i.e. a point where the minimum is attained) of 
r{Q{X)), for r G D[Y], degr < d, on Z{p{Q{X)),K n ), or checking that the 
minimum is not attained and computing the infimum, can be done within the 
same number of operations, and for D = Z within the same bitwise complexity, 
as the computation of Theorem \1.2l 

A proof of the latter, and a number of applications in mathematical program- 
ming, will appear in the continuation of the present paper. 

An easier than optimization problem is the problem of checking whether 
the set Z(p(Q(X)), K n ) is empty, i.e. the feasibility problem. Our immediate 
predecessor here is j2], where it was shown that for homogeneous Q the empti- 
ness of Z{Q 1 (X) 2 + ■■■ + Q k (X) 2 , W n - {0}) can be checked in operations 
in D. 

For the sake of completeness, we state the following straightforward impli- 
cation of Theorem 11.21 

Corollary 1.6. The emptiness of Z(p(Q(X)), K n ) can be checked within the 
same complexity bound as in Theorem U.'A □ 

Remark 1.7. It is easy to see that the bounds of Theorem U.'A are close to 
best possible. Indeed, any real solution of degree 4 equation 

{xl - l) 2 + • ■ ■ + {xl - l) 2 = 0, 

or the system of n quadratic equations 

x\ = x\ = ... = xl = \ 

has coordinates 1 or —1, and there are in total 1 n of them. In the continuation 
of the present paper we will further sharpen this by showing a similar result 
for a system of one cubic and two quadratic equations. 

In a nutshell, the procedure at the core of Theorem 11.21 that we are going 
to describe works as follows. First, we write down the equations for the critical 
points of the projection map X \— > X\ on Z = Z{p{Q{X)),K n ) by equating the 
gradient of p{Q{X)) with the vector proportional to the gradient of X i— > X\, 
that is with a vector of the form (A, 0, . . . , 0). These equations have a rather 
special structure: the variables X occur either within Q{X), or linearly. By 
introducing k new variables (Yi, . . . , Yj,) = Y = Q{X), we thus obtain a system 
of linear equations A{Y)X = b(Y) in X. The next step is to solve this system; 
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we simply loop through all the maximal (by inclusion) candidates for invertible 
submatrices A uw of A(Y) and the corresponding partition X w UX^ of X into 
X\y and the remaining variables X^r. For each of them we rewrite the system 
to express X w as rational functions X w = A'(Y, X^y) of Y and X^- This 
certainly only makes sense, from the complexity point of view, when rk(A(Y)) 
never drops below certain threshold. We make sure by means of an infinitesimal 
deformation that rk(A(Y)) > n — k. Then \W\ < k and the coordinates of X 
are expressed as rational functions in at most 2k variables Y and X-^. We 
are able to describe an isomorphism of a semialgebraic subset, that we call, 
following 2J, piece, of the critical points of X \— > X\ on Z, that corresponds to 
a particular Auw being maximal and invertible, to a semialgebraic subset of 
F fe+ ' H/ ', for F being a real closed extension of K, that is defined by polynomials 
of degree 0(nd). These pieces cover the whole set of the critical points just 
mentioned. 

Finally, we find representatives of connected components of the pieces over 
F, obtaining Y and X^ with values in F, and recover X w and X^ in the 
original field K by computing the limit. 

The actual implementation of this procedure is more involved. Section H] 
describes in detail the candidates for invertible submatrices A uw of A and 
presents the explicit isomorphisms of pieces to semialgebraic subsets in F fc+ ' w ' 
mentioned above. In order to apply the result of Section |U to Z, one needs 
to deform p in such a way that becomes a regular value of p(Q(X)) and of 
p(Y). Further, one needs to deform Q so that the number of pieces of the 
set of critical points of X i— »■ X\ on Z does not exceed (dn) 0( - k \ In fact, our 
deformation will give us a better bound, n°^ k \ on the latter. Lastly, one has 
to ensure (again, using a deformation) that Z is bounded, otherwise we miss 
connected components of Z whose projection on X± is open. 

Our deformations are done by extending K with a number of infinitesimals. 
Subsequent limit computations are needed to recover elements in the original 
set Z by using the following Theorem 11.101 To state it, let us recall some 
notation. For a field F and a transcendental (, we denote by ¥(() C F((£^)) 
the subfield of Puiseux series algebraic over ¥((). For 

(1.8) a = ^a i C 7<? GF((Ci)), 0<qEZ 

with the order v/q > 0, a u ^ 0, define the standard part (cf. e.g. [TU]) of a to be 
a o'i m it is called the limit ao = lim^ a. Note that if v < then lim^ a is not 
defined. When ( is a vector of infinitesimals £i ^> £2 ^> • • • ^> Q, the notation 
lim^a is a shorthand for lim^ 1 (lim ( ^ 2 (. . . (lim^ a) . . . )). It is often helpful to 
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view ( as a parameter and computing lim^ a as computing lim^ a, where lim 
is understood in the usual sense. Note that lim^ is a ring homomorphism of 
the ring ¥(() b = {a G ¥{() | u(a) > 0}, of all the elements of ¥(() bounded 
over F, to F. 

Let F(e) be a real closed extension of a real closed F with infinitesimals 
e = (ei, . . . , Ei) such that £x ^> e 2 3> ■ ■ ■ £t, and let D C F be a computable 
subring of F. For F G D[e][ii, . . . , let Z F = Z(F{Y),¥{e) q ~ 1 ), and let 

^ : F^) 9-1 — > F(£:) m be a rational mapping 



Theorem 1.10. Let F and ^ be as above, with the Y -degree of F at most d 
and the Y -degrees ofVLi and A less than d—1, and their e-degrees at most d. A 
set of univariate representations u(T) of a set of points meeting each connected 
component of lim e ^(Zp) C F m can be computed in (m + d)°^ arithmetic 
operations in D. 

The degrees of the polynomials in u(T) are at most d° iq) . When D = Z, the 
bitsizes of the coefficients of u(T) and of the intermediate data are bounded by 
a polynomial in d, m and d°^ qt> times the bitsize of the input data. 

Theorem I1.1UI generalizes Alg. 11.61] to non-identity mappings 

The remainder of the paper begins with presenting the procedures behind 
Theorem ll.101 along with its proof, in Sections|21and|ni Then Section|Upresents 
the aforementioned decomposition of the zero set of p(Q(X)) under the reg- 
ularity conditions. Finally, Section describes the deformations of p(Q(X)) 
that are needed and completes the proof of the main Theorem 11.21 



As the first part of the proof of Theorem I1.1U1 in this section we address the 
problem of finding limits of the images P(x) of the real roots x G Z(B, ¥(e) q ) of 
a 0-dimensional polynomial system B C ¥(e) q with respect to e — (ex, . . . , ei) — > 
under a polynomial mapping P : ¥(e) q — > ¥(e) m . 

The separating element based methods for finding limits of the real roots 
of B, such as [HI], Chapter 11], cannot deal directly with this situation, 
even in the simplest case of P being the orthogonal projection onto a subset of 
coordinates. 

In this section we generalize these methods to accommodate our needs. 
We introduce P-separating elements a G F(£)[5i, . . . , S m ) such that the map 



(1.9) 




2. Limits of solution images: dimension 
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P(x) I— > a(P(x)) is injective on P(Z), where Z = Z(B,¥(e) q ), that is to say 
that a(P(y)) ^ a(P(x)) whenever P{y) ^ P(x) for x,y G Z. We introduce 
below, in Section 12.2] the corresponding notion for the limit setting, well-P- 
separating elements. 

It turns out that the machinery of ^H], see also pQ, generalizes here rela- 
tively smoothly. Let an ideal (W) C F(e)[S'i, . . . , S n ] = ¥(e)[S] be generated 
by its Grobner basis W C DfejfS 1 ]. That is, we fixed a particular monomial 
ordering on F(e) [S], and the leading (with respect to this ordering) terms of W 
generate the ideal of leading terms of (W). Following Sect. 11.3], we call W 
parametrized special if it is of the form W = {biSf 1 + U\, . . . , b n S^ n + U n } , where 
the leading terms are biSf 1 ,. .., b n S^ n , and deg(£/j) < dj, deg s .(C/j) < for 

1 < i, j < n. Note that the quotient algebra F(e)[5]/(W) has the natural basis 
U(W) of monomials under the staircase, that is, of monomials S a = S^ 1 . . . 
with cii < di for 1 < i < n. In particular, the dimension of the quotient alge- 
bra is d\ . . . d n . In order to keep doing the arithmetic in D[e] when reducing 
with respect to W (and this is one of the purposes of a parametrized special 
Grobner basis), one works in the basis U(W) = {b^S a \ S a e U(W)}, where 

\a\ — aii + h cv„ and 6 a common multiple of &i,. . . ,b n & D[e\. 

The basis W naturally appears when critical points of a coordinate projec- 
tion of a certain special type of hypersurface are computed, as in Sect. 11.6], 
in contrast to a common situation when calculation of a Grobner basis of an 
ideal is computationally very costly (the latter can generally require doubly 
exponential, in the number of variables, running time, cf. [14] ) . Given W, 
one can efficiently compute the multiplication table of the quotient algebra, see 
Alg. 11.22]. Namely, when the degrees of the elements of W in S (resp. 
in e) are bounded by d (resp. by A) it takes (d\) 0( - nt) operations in D, the 
e-degrees never exceed X(nd)°^; when D — Z, the bitsize of the data involved 
is bounded by a polynomial in n and (Xd) n£ times the bitsizes of the elements 
of W, cf. [5, pp. 381-382]. 

Theorem 2.1. For D a computable subring of F, let B C D[e][X 1 , . . . ,X q ) 
define a 0- dimensional polynomial system, that is its own special Grobner basis 
with the LCM of the leading terms equal to bjg G D[e). Let 

N = dim¥(e)[X}/(B) and Z = Z{B, ¥(e) q ). 

Assume the degrees of the elements of B in X as well as in e, and the degree ofbg 
in e, bounded by d. Let P be a polynomial mapping X \— > (Pi(X), . . . , P m (X)), 
with Pi G D[e] [X] of degree < d in X and at most d in e. 
Then a set of at most (m — 1)N 3 candidates for univariate representations 
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u(T) G D m+2 for the elements of lim e P(Z) can be computed in (m + N)°W 
arithmetic operations in D. The degrees of the polynomials in u(T) are at most 



When D — Z ; the bitsizes of the coefficients of u(T) and of the intermediate 
data are bounded by a polynomial in d, m and N e times the bitsize of the input 
data B. 

Note that B will be constructed in Theorem Kill 

Remark 2.2. The number (m — 1)N 3 of candidates can be reduced to N, see 
Remark \2. 29i later in this section. 

The remainder of the section is devoted to the proof of Theorem 12.11 

2.1. On P-separating elements. Here we prepare the ground for the limit 

computations. Denote by x( a ? T) the characteristic polynomial of a linear 
transformation a G A. The Stickelberger's lemma [3 Thm. 4.69] states in par- 
ticular that 



where fi(x) is the multiplicity of root of B. For a given a G A, denote 

by [x] C Z the equivalence class of x with respect to the equivalence relation 
defined by a, so that x is equivalent to y when a(x) = a(y). Then 



Let 7^ b G A be an a- class function on the relation [*] induced by a, that is, 
b(y) = b(x) for any y G [x]. Let S be a variable and consider polynomials 



N. 



(2.3) 



X (a,T) = Y[(T-a{x))M. 



(2.4) 



X (a,T) = ] [ (7 o(i))"W, where /i N = ^ /x(y). 




(2.5) 



X (a + 56, T) = Yl ( T ~ a (%) ~ Sb(x)) 



[x]CZ 



(2.6) 



(a,6,T) = 



Then 



(2.7) <?(a,6,T) = - ^6(x)// N (T-a(x)rw- 1 J] (T-a(j/))"M. 

[x]cz M^Mcz 
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Observe that 

^M-i)(a,6,T)=-6(x)(M [a! ])! ]J (T — a(y)Y [v] + (T — a(x))h(T). 

[*]?%] cz 

for some polynomial h(T). Therefore 

(2.8) g^' l \a ) b ) a{x)) = -b{x){ m )\ \{ (a(x) - a{y))^. 

[*]?%] 

In particular the following holds. 

Lemma 2.9. For a P 'separating a, any O^rG F(e) and x & Z, 

2.10 Pi(x) = — t tt lor each 1 < i < m. 

y ! y ! g^-V(a,r,a{x)) ~ ~ 

Proof. As a is P-separating, gr^ t W~ 1 )(a ) r, a(x)) never vanishes, and [x] ^ [y] 
as soon as Pj(x) 7^ Pi{y)- Thus Pj(X) is a a-class function, and ()2.8|) holds 
for b = Pi. Now (|2.8j) implies the statement of the lemma, as the terms 
r (f i [x]V- l~[[x]^[ y ]( a ( x ) ~ occur in both numerator and denominator of 

the right-hand side of (I2~TU1) . □ 

Note that iV = J2 x ez M^)- To compute the coefficients of x( a + Sb,T), it 
is convenient to write it as 

N 

X(a + Sb, T) = J2 h i ( S ) T ^ where b i e ^[S] , de S b i = N ~ J- 

We need to be able to compute the trace Tr (f(S)) of a linear transformation 
f(S) G A[S]. By additivity of the trace, this is easy to do once a basis U(£>) 
of A, the multiplication table for A in U(*B), that is, a tensor A^ specifying 
linear combinations 

auj = A^i for a, ueX1(B), 
teu(0) 

and the expression of as a linear combination of the elements of U{B) 
with coefficients in ¥[S] are known. Namely, 



Tr(f(S)) = Tr( £ = E E 

wGU(B) aeU(B) ^eU(B) 



a 
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and the computation can be done separately for each coefficient of the polyno- 
mial Tr(f(S)). 

It is well-known (cf. e.g. |5J Thm. 4.69]) that forfeA 

(2.11) Tr(f) = J2v(x)f(xy, j>0. 

xez 

It follows that 
(2.12) 

Tr((a + SbY) = ^2fi(x)(a(x) + Sb(x)y = ^ fi [x] (a(x) + Sb(x)Y , j > 0. 
xez [x]cz 

Then, the bj(Sys (that is, the elementary symmetric functions of the roots) can 
be computed knowing the power symmetric functions (also known as Newton 
sums) Tr(a + Sb), . . . , Tr((a + Sb) N ) of the roots of x( a + Sb, T). Namely, the 
following holds (cf. (4.2), (11.8)]). 

dx(a + Sb,T) _ , qhT ,ST Tr((a + ^n 
— - x{a + ao, i ) Tj^i • 

m>0 

By equating the coefficients of T\ for each j satisfying — 1 < j < N, on the 
both sides of the latter, and recalling that x( a + Sb, T) is monic in T, that is, 
bN^S) = 1, one obtains the following. 

Lemma 2.13. Let x( a + Sb,T) = /^ l j =0 bj(S)T : ' be the characteristic polyno- 
mial of a linear transformation a + Sb £ A[S]. Then 

1 N-i 

(2.14) b i (S) = --—Y t b i+j (S)Tr((a + Sby), 0<i<N-l, b N = 1 

i=i 

gives a recurrence for bi(S) 's, for i = N — 1, N — 2, . . . , 0. □ 

REMARK 2.15. Formulae similar to 112.14)) are known since long time, and 
attributed to \12) . An explicit expression for bi(S) in terms of a determinant of 
certain "almost Toeplitz" matrix with entries specified by Tr((a + SbY) 's can 
be found by using JI3 Ex. 1.2.8]. See also L 13, (2.14)]. 

Using the latter lemma, we can construct x( a , T) and g(a, b, T) given a, b G 
A. We do not need to compute x( a + Sb, T) completely; namely, only ^-linear 
parts of bi(S) and Tr((a + bSY) need to be computed, in view of (|2.14|) and 
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(J2.6j) . To avoid the necessity to handle rational expressions arising from the 
term in (|2.14|) . compute iV! x(a, T) and N\ g(a,b,T) instead. 

In what follows we restrict ourselves to separating elements of the form 
a(P(X)), for a G D[T 1: . . . , T m }. Note that a P-separating a exists and can be 
chosen as follows, for some < j < (m — 1)(^): 

■in 

(2.16) a(P(X)) = a(j,P(X)) = ^f^P^X). 

i=l 

To see this, one proceeds as in e.g. Lemma 4.60]. Let s ^ y G P{Z), and 
observe that the univariate polynomial a(Y, s) — a(Y,y) = YlT=i( s i ~ 1 
is not identically 0, and has at most m — 1 roots. Thus by avoiding at most 
m — 1 values of j, one can make sure that a separates s and y. As there are at 
most (^) distinct pairs of s and y as above, the claim follows. 

Thus we can construct univariate representations of the elements s G P{Z) 
of multiplicity /i(s) = /x + 1 in the form 

(2.17) u(T) = iV! ( X (a, T), ^(a, r, T), ^(a, rP^X), T), . . . , 

^)(a,rP m (X),T)), 

where r G _D[e] is chosen so that the functions ra(X), rP^X),. . . , rP m (X) of A 
are D[£]-linear combinations of the basis elements of A. The latter are chosen 
so that the entries of the multiplication table of A belong to D[e], as dictated 
in turn by the coefficients of the leading monomials in B. Taking r to be the 
LCM bf3 of the coefficients of the leading monomials in B suffices; the numbers 
fi(s) < N are not known a priori, thus we have roughly iV-fold redundancy in 
the output. 

To summarize, we have the following. 

Proposition 2.18. Let B, P be as in Theorem EQ1 A = ¥(e)[X}/(B) be of 
dimension N, and a G Z[T 1; . . . ,T m ] define a P-separating element a(P(X)) 
given by \2. 1 6]) with coefficients in Z of size at most O (log mN). 
Then a set of at most univariate representations u(T) G P[e] m+2 of 

the form A2.17]) . containing for each s G P(Z), a representation \2. 1 7]) with 
fx = /i(s) — 1, can be computed in arithmetic operations in D. The 

degrees of the polynomials in u(T) are at most N, and their coefficients are of 
degree at most 0(d 3 ) in e. 

When D — Z, the bitsize of the coefficients of u(T) and of the intermediate 
data is bounded as stated in Theorem 12. 11 
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The complexity analysis is very similar to algorithms in |3J Chapter 11]. The 
most expensive part is computing the appropriate multiplication table for A, see 
the exposition preceding Theorem 12. II and |5J Alg. 11.22]), and this is identical 
to the special case considered in [loc.cit.]. Note that Pj's are expressed as linear 
combinations of the basis elements of A and therefore our setting for complexity 
analysis of computation of u(T) is essentially the same as in [loc.cit.]. 



2.2. Computation of the limit. We proceed to computing the limits of the 
points given by univariate representations u(T) of Proposition 12.181 We show 
that the limits of points in P(Z) correspond to the limits S <OQ of bounded roots 
of x{ a i T)- Then we normalize the polynomials in u(T) by a Puiseux monomial 
in e that makes the coefficients of x{ a -> T) an d of the rational functions Pi{T) = 
9 ttt^t^P , for each 1 < i < m, bounded. At the same time the values of the 

limit of the normalized denominator T) will be nonzero on S <OQ . Thus 

lim £ Pi(T) at S <OQ can be computed by taking the limits of the coefficients 
of Pi(T), and then evaluating on the elements of S^. We will give explicit 
formulae for lim e Pi(T) in terms of the appropriate repeated derivatives of the 
numerator and of the denominator. 

We need some further notation related to a real closed extension F(£) of a 
real closed field F by infinitesimals £1 3> C2 ^ ■ • ■ d, and its algebraic closure 
F(£). Let O^r 6 F(Q. Then r can be written uniquely as r = £°( r ) (in (t)+t'), 
with 7^ in(r) G F and o(r) G Q e such that ("^ is the biggest, with respect to 
the order in F(£), ^-monomial of r, and r' satisfying lim^r' = 0. In particular 
t' is bounded over F, and lim^ r = in(r) iff o(r) = 0. Further, for 7^ v G F(0 n 
define o(f) = max o(t>,-), with max taken in the sense of the ordering in ¥((). 

l<j<n 

Note that boundedness of r = 3ftr + z3r G F(^) is understood here and 
elsewhere in this section in the usual sense of the norm v/(9fcr) 2 + (3r) 2 being 
bounded over F. As well, lim^ r is understood purely algebraically, that is, 
in the appropriate order setting to the corresponding infinitesimals. In just 
introduced notation, r is bounded if and only if either o(r) = 0, or the rightmost 
nonzero entry of o(r) is positive. 

Let f(T) = Y,7=o c i Ti G F (0P1- Then o(f) is defined to be such that C° (/) 
is minimal, with respect to the order in F(£), monomial making (~ o( -^Cj for 
< j < in bounded over F. In fact o(f) = o(cj) for some j. Define 



(2.19) 



f(T) = \imC° if) f(T), 
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where the limit is taken coefficient-wise. For F(T) G F(£)[T] 9 , we denote 
(2.20) o(F) = max o(Fi). 

l<i<k 

The following statement is an extension of Lemma 11.37 from j3] adjusted to 
the non-multiplicity-free situation, and will be used repeatedly. 

Lemma 2.21. Let f(T) G F(C)[T] be monic. Denoting Z f = Z(/(T),F(C», 
one has 

o(f)= J2 Mr)o(r), Z(f(T),W)=\imZ f . 

rez f , 

unbounded 

Let y G Z(f(T),¥). Then fi(y) = /i(r) equals the multiplicity of y as a 

rez f , 

lim^ T=y 

root of f. Here the summands /i(r) denote multiplicities of roots r of f . 

PROOF. (Sketch.) The coefficients /j of / = SiLo/*-^' wnere fd = 1, are 
elementary symmetric functions of x G Zf taken with multiplicities. Let £ 
denote the multiset of roots of /. Then 



(2-22) = J2 n 

implying 



ecs ree 

|8|=i 



o(fd-i) < max 5^o(r), 

|ef=t -ree 

where the inequality might be struct due to a possible cancellation of higher 
order terms in the sum (|2.22|) . When one of the multisets G equals T, the 
sub-multiset of unbounded roots of /, this inequality turns into equality o = 
o(fd-\y\) = J2 T er°( T )' as ^ ne or der of the remaining summands in ()2.22|) is 
strictly less than o. As o > Xlree °( r ) ^ or an y ® — an d we obtain o(f) = o. 
Therefore 

(2.23) f(T) = lime-° {f) f(T) = lim f[ (T - r) JJ e-° {r) {T - r) = 

rGS-T t£T 

(2.24) =IJ(-in(r)) J] C^-limr), 

t£T T6E-T 

and the first part of the lemma follows. The second part follows from ()2.24j) . □ 
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A P-separating element a will be called well-P -separating (with respect to 
lim e ) if the following conditions hold: 

1. for any s,t/6 P{Z) such that lim e s 7^ lim £ y one has a(lim e s) 7^ a(lim £ y); 

2. o(P(u)) = o(a(P(u))) for any u G Z. 

In particular, |21 implies that if s G P(Z) is unbounded over F then a(s) is also 
unbounded over F. 

Lemma 2.25. Let a,b G D[T±, . . . , T m ] be linear. Let a be well- P-separating. 
Then for any s G \rm £ P(Z) of multiplicity //[ s ] 

g^-^q^ais)) g^^ja, P 3 ,a(s)) 
b(s) = —, rr- — -. In particular, s,- = —7 tt- — — — , 

for any 1 < j < m. 

Proof. As a is well- P-separating, b is an a-class function on P(Z) as well as 
on lim e P(Z). Note that b satisfies lim e b(P(x)) = 6(lim £ P(x)), for any x £ Z 
for which lim e P(x) is defined. We shall express lim e 6(P(X)) as a univariate 
rational function, that gives lim £ b(P(y)) when evaluated at a(P(y)). With 
a = a(P(X)) and b = b(P(X)), denote 

r](a, b, S, T) = lime- o(x{a ' T)) x{a + Sb, T). 

e 

Denoting by Z <00 the set of x G Z such that a(P(x)) is bounded, using (|2.5j) 
and Lemma [2.211 fin particular (|2.24j) ) one obtains 

n{a, b, S, T) = \ime-° ixia ' T)) JJ (T - a(P(x)) - S6(P(x))f w = 

[i]Cl 

= Y[ lim(T-a(P(x))-Sb(P(x))ywx 

[x]CZ <00 



] [ lim (e- o(a(p(x))) (T - a(P(x)) - Sb(P(x))))' 



_ _ e 

[x]CZ-Z <x 



= G(S) Yl (T -a(y) - Sb(y))^, where 

[y] Clime PfZ) 

= |"J (-in(a(P(z))) - S&(lini£-°( p ^P(a;))f w G F[S], 

[x]CZ-Z<oo 
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as o(a(P(x))) = o(P(x)). Moreover, G(0) ^ 0, as m(a(P(x))) ^ on un- 
bounded x G P{Z). 

In view of (|2.24|) one obtains in particular 

(2.26) x(a,T) = C J] (T - a(y))^, with 

[y]C]im E P(Z) 

where denotes the multiplicity of a(y) as a root of x(a, T). 

Next, we compute, for g(a,b,T) defined by (|2.6|) . lim e £~°( x(a ' T ))(yf(a, 6, T). 
We see that it equals to g(a, b, T) (if it would not be the case, it had to vanish 
identically, as o(x(a, T)) > o(g(a,b,T))), as defined in (j2.19j) . so we get 

g(a, b, T) = JLr,{a, 6, S, T)\ s=0 = -G'(0) J] (T — a(y))"M- 

[y]C\imeP(E) 

-G(0) £ ^^(^(T-a^))^" 1 J] (r-a(s)rw. 

[j/] Clim e P(Z) [j/] ^ [s] Clim e P(Z) 

Then, for any y G lim e P(Z), 

$<"M-i>( fl , 6, T) = -G(0)&(y)(/i M )! J] ( T " a ( s )) MW + ( T " %))#( T )> 

[s] Clime P(Z) 

for #(T) G F[T]. Hence 

^w-^Co, 6, o(y)) = -G(0)6(y)(p M )! J] (o(y) - o(a))"w 

[y]^[s]Clim £ P(Z) 

and we obtain, in view of (|2.26|) . the statement of the lemma. □ 

Let us show that a can be taken to be a = a(j, P(X)) for a certain j as in 
(|2.16j) . The only difference with the argument above is that we have to avoid 
more "wrong" values of j. Let x, y G P(Z) be such that lim £ x and \im £ y 
exist and are not equal. Then the polynomials a(Y,\im £ x) — a(Y, lim £ y) and 
a(Y,x) — a(Y,y) are not identically and each of them has at most m — 1 
roots. Thus by avoiding at most 2(m — 1) values of j, one can make sure that 
a separates lim e x and lim e y, as well as x and y. 

To ensure the remaining condition in the definition of well-P-separating 
element, consider W = {hm £ e~°^s ^ s G P(Z)}. Choose j such that 
a (j, w ) f° r an y w W. Such a choice is always possible: a(Y,w) G ¥[Y] 
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has at most m — 1 roots; thus avoiding (m — values of j achieves the 

required. Then o(a(s)) = o(s) for all s € P{Z), implying condition |2] of the 
definition. 

As there are at most (^) distinct pairs of s and y as above, and since 
\W\ < N, we obtain 

Lemma 2.27. There exists an integer < j < (m-l)iV 2 such that a(P(X)) = 
a(j,P(X)) is well- P -separating. □ 

Combining Lemmas I2.25l and l2.27l gives for an appropriate a at most 0(\/N~) 
candidates for univariate representations u(T) for the points in lim E P(Z), as 
x(a,T) has at most 0(y/N) different root multiplicities. We outline now how 
u(T) are actually computed. Let b = Pi for some 1 < i < m. We loop 
through 1 < j < (m — 1)N 2 in order to be sure to find an appropriate well-P- 
separating a(P(X)) = a(j,P(X)). This means that we will return candidate 
representations for each such j. 

A further technical point is that we operate in the ring D[e] rather than in 
a field. We utilize the idea of Remark 11.44] for r = b^: 

(2.28) x(ra + Srb,rT)=r N x(a + Sb,T) 

and proceed similarly to the procedure of Alg. 11.45]. 

Remark 2.29. Compared to (loc.citj, our simplification is that we do not 
try to filter out "wrong" a; such a check would require finding \P(Z)\ to be 
able to verify that a is P-separating, and then proceeding similarly to remarks 
on p. 398]. To obtain "good" a, one would first select a's with the biggest 
degree of x{a, T); among the latter select a's that are P-separating, by choosing 
a's with minimal degree of gcd(x(a, T) , dx ^j ) ■ To ensure that a is well-P- 
separating, one would select a's with minimal degree of gcd(%(a, T), d ^ T ^ )- 
Finally, to make sure that o(a(x)) = o(x) for any x G P{Z), one would check 
that o(%(a, T)) > o(g(a, Pi,T)) for all 1 < i < m. 

Anyhow, this shortcut does not worsen the asymptotic running time. 

Thus we compute x( ra + Srb,T), as already outlined in the first part of 
this section, and then operate with e~ o( - x ^ ra ' rT ^x( ra + Srb, rT) to obtain the 
remaining data for the limit computation. We have 

X (ra + Srb, rT) = lim£-° (x(m ' rT)) x(ra + Srb, rT). 

£ 

Using the latter to compute g(ra,rb,bT) and g(ra,r,bT) as in the proof of 
Lemma EM by (|2~28jl we obtain b(x) = ffej^gff for each x 6 lim e P(Z) 
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of multiplicity fi + 1, as required. For each j we loop through all the possible 
values of //. Thus in total we will have no more than (m — l)iV 3 univariate 
representations, as stated in the theorem. 

The complexity analysis for a given j runs parallel to the analysis given 
in Sect. 11.5] (see also a remark following Proposition 12.181 above) and is 
omitted. This completes the proof of Theorem 12.11 

3. Limits of solution images: dimension > 

In Section |2] we described a procedure to compute limits of images of finite 
algebraic sets under polynomial mappings. To complete the proof of Theo- 
rem ^^]J we proceed to computing limits, with respect to e — > 0, of rational 
images of samples of connected components of real algebraic sets, by reducing 
to the 0-dimensional case. 

For F G D[s][Yx, . . . let Z = Z(F (Y),¥(e) q - 1 ), and let * be as in 

(jl.9|) . We want to find points in each connected component of lim e ^(Z ) C F m . 

Theorem 3.1. Let Fq and ^ be as above, with the Y -degree of Fo at most 
d and the Y -degrees of Qi and A less than d — 1, and their e-degrees at most 
d. Then one can construct a (q + 2)-variate O-dimensional polynomial system 
B over the real closed transcendental extension F(ex, . . . ,££+2) of¥(e) by two 
infinitesimals < Ei +2 £e+i £t, with zero set Z(B) such that 

C = lim #(tt(Z(B))), 

ei,...,£f + 2 

where ti is the orthogonal projection to Y±, . . . , Y q _i, intersects each connected 
component of lim eit Ee ^/(Z ). 

The system B can be constructed in (m + d)°^ operations over D, and it 
satisfies the properties required by Theorem \2.1l the set C can be computed 
using the procedure of Theorem \2.1l 

In the proof we will denote \i = ee + i and ( = ee+2- 

First of all, let us reduce the setting to the case when A = 1. Assume that 
Fq(Y) > 0, otherwise replace F with Fq. Consider 

F(Y l7 ...,Y q ) = F (Y) + (1 - Y q A(Y)) 2 , Z = Z(F(Y),F(e)*). 

Then V(Z ) = P(Z), where P G D[e]\Y u Y q ] m is given by 

(Y 1 ,...,Y q )^(Y q Q 1 (Y),...,Y q n m (Y)). 
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Thus from now on we consider the problem of computing representatives of 
connected components of the set lim e P(Z) C W m . Note that the map P is as 
required by Theorem 12.11 

The idea is to look for points of minimal norm in the connected components 
of P{Z). We will repeatedly use the following statement. 

Proposition 3.2. ( L 5, Prop. 11.56]) Let S C K(e) u be a semialgebraic set, 
for K(e) a real closed extension of a real closed field IK by an infinitesimal e. 
Then lim £ S C IK n is a closed semialgebraic set. If S is in addition connected 
and bounded over IK then lim e S is connected. □ 

Introduce another variable Yq that will help us in this task; we will write 
down, after necessary preparations, equations that specify the critical points of 
the projection Y i— > Yq. Replace F by 

m 

F'(Y , ...,Y q )= F(Y) + (Y -J2 mff- 

i=i 

We extend P to the zero set Z' of F' by "ignoring" Yq. Obviously P(Z) = 
P(Z'). 

Next, we need to make the set Z' bounded, in the standard way of e.g. jSJ 
Chapt. 11]. Introduce an infinitesimal < C e<, a variable Y q+ \ and define 

9+1 

P,(Y , • • • , Y q+1 ) = F'iY) + (1 - /i 2 Yf fi %T = Z ^ ¥ ^ V) q+2 )- 

3=0 

Again, extend P onto Z<°° by "ignoring" Y q+ i. Define 

Z, = Z(F',¥(e,^ +1 ) 



Lemma 3.3. The following holds: 

(i) P(Z<°°) is closed; 

(ii) P(Z<°°) C P(Z,); 

(in) lim At P(Z<°°) equals the closure of P(Z). 
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Proof. The first part of the lemma follows from [5J Thm. 3.20], as P(Z<°°) 
is the image of a closed and bounded semialgebraic set under a continuous 
semialgebraic function. 



The second part is straightforward, by observing that the projection of Z 



<oo 



on the first q + 1 coordinates is a subset of Z^. 

Then, lim M P(Z< 00 ) is closed by Proposition E21 Next, lim M P(Z<°°) D 
P(Z). Indeed, if u = P(y) with y E Z then there exists G F(e,/i) such 
that /i 2 |b/J 2 = 1, where we denoted y^ = (y Q , y x , . . . , y q+1 ) G Z<°°, implying 
FM = and it = P{y„) G P(Z<°°). 

Finally, the inclusion of lim^ P(Z<°°) in the closure of P(Z) follows from 
the second part. □ 

The next step is to introduce a deformation, same as in Sect. 11.6], 
that ensures smoothness. As in jJJ (11.13)-(11.14)], introduce an infinitesimal 

< ( <C /i, and set 

(3.4) F ( (Y) = c(^ 1 (Y; d + ^(Yf + Y^ - 2q - 3^ + (1 - QF»{Y). 

Note that the difference between |3J (11.13)-(11.14)] and (|3.4|) is purely nota- 
tional: variables and infinitesimals have different names, and they range slightly 
differently (i.e. there are q + 2 variables Y instead of q variables X). Here d 
is the minimal positive even number strictly bigger than the degree of F^; this 
simplification of p, (11.13)-(11.14)] is explained in Rem. 11.49]. 

By 11.50-11.51], the set Z c = Z(F C , F(e, //, C) 9+2 ) is a nonsingular alge- 
braic hypersurface contained in the ball of radius -, and such that lim^ Z$ = 
Z<°°. In particular lim^ u exists for any u G Z^. 

Consider the set of the critical points of the projection Y \— > Y on Z^. 
They satisfy the system of equations 

(3.5) = P C (F) m m 



dYi dY Q 



The zero set ~D of ()3.5|) in the algebraic closure of F(e, /i, C) is analysed 
in detail in Sect. 11.6]. In particular, \K ( \ < (cf. [5, Prop. 11.57]). 

Note that F^(Y) can be reduced (cf. [51 Notation 11.59]) using the relations 

OF — 

gy- — 0, j > 1 so that the degree in Yj of the result is strictly less than d. The 

^ OF, 

result of these reductions, together with the j > 1, forms the polynomial 
system B of Theorem 13.11 we are proving. 
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To complete the proof that C is as claimed, it remains to prove the follow- 
ing lemma. It will make sure that we recover in £ a representative of each 
component of lim e P(Z). 

Lemma 3.6. The set lim^^ P(K^) intersects each connected component of 
P(Z) in a point of minimal norm. 

PROOF. By construction, P{Kq) contains the local minima of Y i— > Yq on 
P(Z<). 

By 11.55], the set = lim^ (the set of pseudo-critical points) inter- 
sects each connected component of Z<°°. 

As P does not depend upon (, we have that 

limP(Z ? ) = P(limZ c ) = P(Z<°°). 

Due to closedness and boundedness over ¥(e,fj) of P(Z<°°) (cf. Lemma l3~3l . 
Y i — ► Yq reaches its minimum on each connected component C P(Z<°°). 
Note that by Proposition 13.21 we have = lim^C^, for a connected Cq C 
P(Z^). Consider P _1 (C /1 ) C Z<°°. This set is a union of connected components 
of Zp°°, and intersects each of them; among these intersection points there 
are ones with minimum value of Yq. Therefore P(Kfj) contains representatives 
of the local minima (that need not be singletons any more) of Y i— > Yo on 
P(Z»). 

As P(Z<°°) is the intersection of P{Zy) with P(Pi/ M ), where By^ denotes 
the ball of radius - and centre at the origin, each connected component C 
of P{Z^) that intersects P^Bu^) will contain a connected component C M of 
P(Z<°°). 

We assume now, without loss of generality, that C contains a point of F(e)- 
finite norm. Consider the intersection C r of C and the (smaller) set P(B r ), 
where 1 + mm ueC \\u\\ < r £ F(e). This is a bounded over F(e) semialgebraic 
set, with finitely many connected components C\. By Proposition 13.21 the set 
lim^C* is connected, for any i. As C % r C C^, and as Yq has a local minimum 
on C*, each C\ fl P{Ky) ^ 0. Thus lim M P(i^) intersects each connected 
component of P(Z) in a point of minimal norm, as required. □ 

A straightforward count of the number of operations needed to construct B 
completes the proof of Theorem 13.11 

To complete the proof of Theorem 11.101 a Pphy Theorem 13.11 (with £ = 2) 
to Zq := Zp and F Q := F. We obtain a 0-dimensional system B of equations 
in the real closed extension of F(e) by two extra infinitesimals jj, — £3 and 
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( = £4, with the set of solutions such that lim^ ^(ir(K^)) intersects each 
connected component of ^(Z Q ) in a point of minimal norm. Hence applying to 
B the procedure of Theorem l2.1l will produce a set lim e A1 ^ ty(ir(K^)) intersecting 
each connected component of lim e ty(Zp) in required running time, as iV < d q . 

□ 

4. Pieces of extrema of a bounded level set over a 

quadratic map 

Let the quadratic map Q : K n — > K k be given by 

Qj :X^ l -X T H 3 X + b^X + c v l<j<k 

Cj e D, b 3 e D n , Hj = Hj e D nxn . 

Let ( G K. We consider the (-level set V = Z(p(Q(X)) - £ K n ) of a 
polynomial p G -D[Y"i, . . . , Y^] of degree d over a quadratic map Q : K n — ► K fe 
given by (|4.1j) . 

We assume that D is a computable subring of a real closed field K and V is 
bounded over K. The latter assumption is in fact technical. It can (and will) 
be ensured in Section |21 by introducing an extra infinitesimal £0 an d two extra 
variables, to make V the image of an algebraic set on the sphere of radius 1/eq 
under projection onto the first n coordinates, and subsequently removing e by 
restricted elimination. 

As well, we assume that ( is a regular value of p(Q(X)) and of p{Y). 

Let V c C V denote the set of critical points of the projection map X (—> Xi 
from V to K. Then, as ( is not a critical value of p(Q(X)), V c is an algebraic 
set defined by 

(4.2) P(Q(X)) - C = 0, 

dp(Q(X)) dp(Q(X)) 
1 ' ) dX 2 dX n 

Due to the assumption that V is bounded, V c intersects nontrivially each con- 
nected component of V, see [TU], Prop. 7.6]. Thus a set S c of representatives 
of each connected component of V c will also intersect each component of V. A 
useful property of S c is that it will contain points in V with minimal value of 

Define 

Pi (X) = ^-, for l<j<k. 
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The set of m x n matrices Ai, . . . ,A k over K, m < n, is said to be in r- general 
position with respect to if 

fe 

(4.4) vk(J2tiAi)>r for all t G (Pi(Q(V)), ■ ■ ■ ,Pk(Q(V))) C K fc . 
i=i 

Note that t is never here, as ( is not a critical value of p(Y). This is a weaker 

notion of general position than the one used in , where t was allowed to range 
k 

over K , instead of belonging to the image of V under the polynomial mapping 
X^( Pl (Q(X)),...,p k (Q(X))). 

Note that all the nonzero K-linear combinations of almost any k-set ofmxn, 
in < n, matrices over IK have rank at least m + 1 — k. 

Below we assume that the matrices Hj's, that are obtained from H/s by 
removing the first row, are in r-general position with respect to V. For U C 
{1, . . . , m} we denote U = {1, . . . , m} — U. For U C {1, . . . , m} and W C 
{1, . . . , n) we denote by A uw the submatrix of A obtained by removing all the 
rows in U and all the columns in W. 

To compute S c in less than exponential in n time, or even just to give an 
upper bound on its size of order less than exponential in n, standard methods 
such as one in [TU], jSJ Chapter 11] that treat the system ()4.2|) - (j4.3J) directly do 
not suffice. 

Indeed, in these methods the number of variables, n in this case, appears 
unavoidably in the exponent of the bounds. In contrast, we are able to get 
better results for r being close to n. We cover V c by (at most n°( n ~ r >) semial- 
gebraic sets we call pieces, each of them isomorphic to a semialgebraic set lying 
in K* with t < k + n — r. (Here covering means simply that the union of the 
pieces is V c ; they in general intersect, and can even be equal one to another.) 
Each of the latter is defined by at most 0(dn 2 ) polynomials of degree at most 
0{nd). To them, one can apply the standard technique, see e.g. to bound 
the number of its connected components, so there will be at most (nd) 0( - k+n ~ r ^ 
of them, and to find their (perhaps non-unique) representatives. However, for 
the latter, for a technical reason, namely the necessity to take limits with re- 
spect to certain infinitesimals (in particular ( will be treated as such), we shall 
use the approach presented in Sections 121 and El 

The main results of this section are as follows. 

Theorem 4.5. Let ( be a regular value ofp(Q(X)) andp(Y), and the level set 
V = Z(p(Q(X)) — (, K n ) be bounded over K. Further, let the matrices Hj 's be 
in r-general position with respect to V . Then one can construct a covering of 
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the set of critical points V c CV of the projection map X i— > X\ by semialgebraic 
sets V C (U, W) indexed by invertible submatrices of J2jPj(Q(X))Hj with row 
sets U C {2, . . . , n} and column sets W C {1, . . . , n} that satisfy r < \U\ — 
\W\ < n — 1. For each such W, the polynomial mapping 



( 4 - 6 ) Y (Q{X) 



w : K n -> K k+n ~ lwl 
( X W 

on each <f> w {V c (U, W)) has explicitly given inverse 

<Puw ■ K k+n - lwl -> K n 

Y\ (Q uw ( Pl (Y),..., Pk (Y),T) 

t) ^ V T 

wiere 0[/vk js a vector of rational functions in Pj(Y) 's and T, all with the same 
denominator. The degrees of the latter and of the numerators are at most \ W\. 
The set (j)(V c (U, W)) is defined by p(Y) = ( and (n - \W\) 2 + k + l + n-r 
polynomial equations and one un-equation in the Pj 's, Y, and T, of degree 
0(\W\). In total there are at most n 0<n_r ) pieces. 



Remark 4.7. The statement of Theorem 14. 51 holds un-amended in the more 
general setting of p being a differentiable function, and "semialgebraic set" 
replaced by "joint zeros of equations and non-zeros of inequations" . The proof 
we give goes through in this case, too. We chose to remain in semialgebraic 
setting for the sake of clarity only. 

The following summarizes, in the setting of the paper, the necessary com- 
plexity estimates for the procedure outlined in the course of proving Theo- 
rem ESI 

Proposition 4.8. In notation of Theorem \4.5\ let K = h(e), with e = 
(ex, ...,££) infinitesimals, £ > 1, D = D^[e], and p G D[Y] of degree d in 
Y. Let the degrees of Q and p in e be at most d! . The map 0^ and the 
definition for (p(V c (U, W)) can be computed within (n(d+ d')) ^ +k+n ^ r "> arith- 
metic operations in D^. In the case = Z, the bitsizes of them and the 
intermediate data are bounded by the bitsize of Q and p times a polynomial in 
n, k, logd, log(l + d'), and £. 
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In the course of the proof of Theorem 14.51 we construct the explicit maps 
4>jjw an d semialgebraic definitions for V C (U, W). Equations (|4.3|) can be written 

as 

k 

dx (x) = H^ x )^^ x + h i) = °> 2 ^ 

where is the standard i-th coordinate vector. Thus, denoting 

k k 

i=i i=i 

where 6, denotes 6j with the first coordinate removed, one can write ()4.3|) in 
the matrix form, as follows. 

(4.9) HQ(X))X = b(Q(X)). 

At this point we would like to give an outline of the remainder of the proof of 
Theorem 14.51 We will compute a set of solutions of the system of equations 
()4.2|) - (j4.9j) that intersects each connected component of V. Doing this in the 
standard way would take an exponential in n number of operations, and this 
is exactly what we want to avoid. 

The structure of the equations (|4.2j) - (j4.9|) suggests the substitution Y = 
Q(X). It turns ()4.9|) into a system of linear equations $>(Q(Y))X = b(Y) 
with respect to X. We cannot simply "invert" $(Q(Y)), as it need not be 
of full rank. However, rk($(Q(y))) will always be at least r, allowing us to 
"split" the solving into fi°' n " r ' cases, one for each maximal invertible submatrix 
(parametrized by U and W), that will be inverted. This gives (for each such 
case) rational expressions 4>iw f° r r °f ^ ne "^' s ^ n terms of Y and the remaining 
n — r of X's (that we will denote T), as well as giving extra (in)equations 
involving Y and T. The latter define the sets V C (U, W), essentially completing 
the proof. 

We proceed with the detailed proof now, preparing the ground for intro- 
ducting in (J4.18|) the variables Y. The r-general position assumption implies 

rk($(Q(x))) > r for any x G V. 

Thus there are at most n°' n ~ r ' maximal by inclusion invertible submatrices of 
$(Q(X)). Indeed, there are at most s(r) = (J2m=r Q)) 2 °^ them, by counting 
number of pairs (U, W) of subsets U C {2, . . . , n}, W C {1, . . . , n} satisfying 
r < \U\ — \W\. If r < n/2 then s(r) < 2° (n) < n 0( - n ~ r \ Otherwise one has 
s(r) < (n — r ) 2 ( n ™ r ) 2 < n°( n ~ r \ again as required. 
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As well, at least one <&{Q(X)) UW will be invertible. Hence 

det ${Q{X)) V , W > = for all \U'\ = \W'\ = \U\ + 1 
(4.10) U C U' C {2, . . . , n}, W r cW'c{l,... ) n} 1 

det$(Q(X)W^0. 



Noting that (ETTUJ) implies that det $(Q(X)) a w 
bigger than \U\, one obtains the following. 



for all U' and W of size 



Lemma 4.11. 
(4.12) 



Vc. 



u 



WW, 



UC{2,...,n} 
W<Z{l,...,n} 
r<\U\=\W\<n~l 



where V C (U, W) is defined by the equations \4.2\) - A4.3\) and by &4.1U\) . The 
number of elements in the union in A4.12\) is at most n 0(jl ~ r \ 

Without loss in generality U = {2, . . . , r + 1}, W = {1, . . . , r}. Invert the 
matrix \1/ = §(Q(X))uw using the Cramer rule: 



(4.13) 



det 
det# 



for 1 < i, j < r, 



where *f>(i,j) is the matrix obtained from $ by removing z-th row and j-th 
column. Then the system <&(Q(X))X = b(Q(X)) can be rewritten in the block 
form as 



(4.14) 



x w 




V 






b w _ 



where the common "(Q(X))" is dropped for the sake of readability. Applying 
( ^-i ?) to both sides of (|4.14|h one obtains 



(4.15) 











u,w 



x w 









$77 w*" 1 ^ 



Thus the following, together with (|4.10j) and ()4.2|) . provides another definition 
of V C (U,W). 

(4.16) X w = HQ(X))^ W ■ (b(Q(X))u - <S>(Q(X)) uW X w ) , 

(4.17) b(Q(X)) w = $(Q(A) )jj w ^(Q(X))^ w • b(Q(X)) v . 
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Observe that in the latter definition of V C (U, W) the only place where X w 
appears other than as an argument to Q is the left-hand side of ()4.16|) . Set up 
the mapping <Pu\y as follows. 

( 4 - 18 ) m , . ■ (b(Y)u ~ HY) uW T) 

One establishes that 4>ijw ac ^ s as claimed in the statement of the theorem. 

Lemma 4.19. The mapping (fiw restricted onto V C (U, W) has inverse 4>uw> that 
is, 

<Puw(MVc(U,W))) = V c (U,W). 

Proof. We have to check that <j>uly(<t>w( x )) — x f° r an Y x e Vc(U, W). As 
neither (f>w nor (f>^ w change x^r, in view of (j4.18|) it suffices to check that 
$>(Q(x))u W (b(Q(x))u — $>(Q(x))uw x w) = x w- But the latter holds as it is 
nothing but ([4.16)1 . a part of a definition for V C (U, W). □ 

We proceed to write down an explicit definition for <p w (V c (U, W)) in terms 
of variables Y and T used in (14. 18)1 . Denoting Q = det$(Y) uw , one has the 
following polynomial equations 

(4.20) p(Y) = C 

(4.21) Q 2 Y = Q 2 Q(<f ) ^ w (Y,T)) 

(4.22) Q b{Y)jj = ^(Y) Vw ^{Y% l w ■ b(Y)u, 

where puzzlingly looking multiplication of both sides of (|4.21jl and ()4.22j) by Q 
clears denominators coming from (|4. lHjl in the right-hand sides, and 

det $(Y) UIW , = for all \U'\ = \W'\ = \U\ + 1 

(4.23) U C U' C {2, . . . , n}, W C W C {1, . . . , n}, 

det$(y)^^o. 

Apart from (|4.21j) . that "bootstraps" Q, these (in)equations already appeared 
above, with K substituted for Q(X) and T substituted for X-^f. 

Lemma 4.24. The relations A4.2U\) - A4.2!^) provide a semialgebraic definition of 
<for(V c {U,W)). 



Computing over quadratic maps I 27 



Proof. Let (y,t) belong to the semialgebraic set defined by ()4.2U|) - (j4.23j) . 
and x = (j)vw{y,t). We shall check that x G V C (U,W). Due to (|OT|) . the 
equation (|4.2|) holds for X = x. Similarly, the remaining (in)equations (|4.1(J|) . 
dnBD - gTTD defining V C (U, W) hold for X = x. 

By inspection, any x G V C (U,W) gives rise to (y,t) = (pw( x ) m the set 
defined by □ 

The entries of the matrix are linear polynomials in p±, . . . ,Pk- Thus 

the determinants of its m x m-submatrices, that come via (|4.13j) into the def- 
inition of V C (U, W) by ()4.20|) - ()4.23|) . will be polynomials of degree at most m. 
Similar straightforward degree counts complete the proof of Theorem 14.51 

We proceed to prove Proposition 14.81 The only nontrivial part concerns the 
complexity of computing the map _1 (£7, W) using the Cramer rule (j4.13|) . 

To count the number of arithmetic operations in required to compute the 
determinants, one can either slightly extend [21 (2.8)], or refer to [3 Alg. 8.38, 
Rem. 8.39(b)], to obtain that computing the determinant of a submatrix of 
$(Y) can be done in arithmetic operations in Dl[e, Y], and then refer to 
[5, Alg. 8.38]. 

Using the latter source, one sees that for = Z, the bitsizes in the answer 
and in the intermediate data will be bounded by (r + log n)n + ( k + 1) log(n(o? + 
d') + 1), with t a bound on the bitsize of coefficients in the entries of &(Y). 
Noting that r is bounded by log d times the bitsize of p and Q completes the 
computation of the bitsize bound for determinants. The remainder of the proof 
is a straightforward use of the complexity analysis of arithmetic operations in 
polynomial rings in [HJ Algs. 8.8, 8.10]. 

5. Proof of Theorem 11.21 

Let Eq > be an infinitesimal over K. We use it to deform Z = Z(p(Q(X)), K n ) 
so that it becomes bounded. (A priori this deformation is not necessary if Z 
is known to be bounded in the first place.) We can assume that p(Y) > for 
all Y, otherwise we can replace p by p 2 . Introduce extra variables Xq and Yq, 
and set 

p(Y)=Y 2 +p(Y). 

Set 

n 

Q (X) = l-e 2 J2 x " 

i=0 

and abuse slightly the notation by setting Q = (Q , Qi, . . . , Q k ). Further, set 
k = k + 1 and n = n + 1. Then for Z = Z(p(Q(X)),K(e ) n ) one sees, by 
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using [SJ Prop. 11.47], that the projection to the last n coordinates of any set S 
meeting every connected component of Z meets every connected component of 
Z(p(Q(X)), K(e ) n ). Assuming one can compute S as a set of univariate repre- 
sentations in D[e ,T], one then can use restricted elimination 5, Alg. 12.43] to 
replace Eq by a sufficiently small element of the field of fractions of D to obtain 
univariate representations of points of Z. 

Next, we deform Q by defining Q(t) : K n — > K[t] k as follows 

(5.1) Q 3 (t, X) = Qj{X) + ^X T diag(P\ 2 j , h j )X, < j < k. 

Obviously Q(t) defines, as well, a quadratic map F n — > ¥ k for any field F D 
K(t). The following lemma states that the Hessians of the Qj(t)'s are in gen- 
eral position, in sense that their nonzero linear combinations (that will be the 
matrices A(Y) mentioned in the discussion in the beginning of this section) are 
of maximal possible rank. The statement is similar to (3.6)]. 

Lemma 5.2. Let the matrix A(Y,T) with the entries in K[Y,T] be defined by 

k 

A(Y, T) = J2 YAHj + T diag(P- 1 , V~\ n^ 1 )). 

5=1 

Let t be transcendental over IK. Then for any held F D K(t) and ^ y G F fc , 
the rank of A(y, t) is at least n — k + 1. 

There exists < t' G IK such that for any i G K satisfying < i < t' and 
the rank of A(y, l) is at least n — k + 1. 

Proof. Consider 

k 

B = B(Y, T, fi) = ]T Y 3 ^H 3 + Tdiag(P- 1 , V~\ n^ 1 )) 

and the homogeneous, with respect to Y, as well as with respect to {/i, T}, 
ideal J = (detB uw | U, W C {1, ...,n}, \U\ = \W\ = n - k}) in the ring 
D [Y, T, fj] C K[y, T, fx) . 

Note that every (y*, £*) 7^ satisfying rk(A(y*, t*)) < n — k + 1 gives rise to 
^ (y*,r, 1) G Z(J) = Z(J,K fc+2 ). Vice versa, ^ (y*,r, 1) G Z(J) obviously 
implies rk(A(y*, z*)) < n - k + 1. 

The idea is now to show that there exists a nonzero polynomial f(T) G K[T] 
such that /(t) = for all (y*,t) ^ with rk(A(y*,t)) < n - k + 1. As t is 
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transcendental over K, that is, it cannot be a root of a polynomial in K[T], this 
will imply the statement of the lemma. 

The ideal J' = (J : (Y)°°) PlK[T, p], obtained by "projectively" eliminating 
Y from J, is homogeneous. By the "Main Theorem of Elimination Theory" , 
see e.g. jTJ Theorem 14.1], the image of Z(J) under the corresponding pro- 

2 

jection is Zariski-closed. Hence it coincides with Z(J') = Z(J',K. ). More- 
over, J' is nonempty, as Z(J), and hence Z(J'), do not contain elements with 
p = and T = 1, as follows immediately from the properties of the diago- 
nal matrices diag(F, 2 J , . . . , Thus J' contains a homogeneous polynomial 
f(T, p) G K[T, fj] that is not divisible by p. 

Hence any t for which there exists y* ^ satisfying rk(A(y* , t)) < n — k+1, 
satisfies f(t, 1) = 0, that is, t is algebraic over K, a contradiction showing the 
first part of the lemma. To obtain the second part, observe that f(T, 1) vanishes 
only on finitely many elements of K. Choose t' to be the closest to root of 
f(T, 1) among the positive elements of K, if such a root exists. Otherwise 
choose i' = 1. □ 

Next, we shall perturb p(Q(X)) so that is not its (and neither that of 
p(Y)) critical value by subtracting an appropriate constant r from it. (Such 
t is called a regular value of p(Q(X)) and of p(Y)). We will talk about the 
r-level set of p(Q(X)), that is just Z(p(Q(X)) — r, K n ). The following is an 
immediate consequence of the semialgebraic Sard's theorem Thm. 9.6.2], [SJ 
Thm. 5.57]. 

Lemma 5.3. Let ¥ be a real closed field and r a transcendental infinitesimal 
over F (respectively, r > a sufficiently close to element of Fj. Then r (and 
any i satisfying 0<i<r)isa regular value of any nonzero f(Y) G ¥[Y]. fn 
particular, provided that ¥ contains the field generated by the coefficients of p 
and Q, one has that t (and any i as above) is a regular value of p(Q(X)) and 
ofp(Y). □ 

Let Eo S\ ^> E2 > be two more extra infinitesimals over K, and denote 
Q = Q(e 2 ). We deform Z as follows: 

Z = Z(p(Q{X))-e 1 ,K(e ,e 1 ,e 2 ) ii ). 

At this point we are ready to use the tool from Section where it is de- 
scribed in slightly greater generality. According to Theorem 14.51 we have a 
covering of the set V c of the critical points of X t— > X on Z by semial- 
gebraic sets V C (U, W). Moreover, Theorem 14.51 gives us for each V C (U, W) an 
isomorphism <fiy/ (given by polynomials in D[X] of degree at most 2d) so that 
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4>w(V c {U,W)) C K{e ,e 1 ,e 2 ) o{k) , as well as its inverse <fiuw (gi ven by rational 
functions, with common denominator, with coefficients in D[e], of degrees at 
most h). By Proposition 14.81 this data can be computed by (n(d + d'))°^ 
arithmetic operations in D. 

The sets 4>w(Vc(U, W)) and V C (U, W) are both defined by equations and one 
inequation A ^ 0, with A G D[e][Y] (respectively, A G D[e][X]). By adding 
one extra variable as in the beginning of Section |3] we convert each of them 
into a real algebraic set: add equation Y^A = 1 (respectively, Ajj +1 A = 1) 
and extend the maps <pw and 4>uly by "ignoring" these extra variables. Apply 
to Zp := (f>w(V c (U, W)) and ^ := the procedure of Theorem II. 101 It will 
produce a set R(U, W) of univariate representations of points intersecting each 
connected component of lim £ V C (U, W). 

By the following lemma, the union of the R(U,Wys over U, W will in- 
tersect each connected component C of Z, as by Proposition 13.21 one has 
C = lim £l . e2 C e , where C £ is a connected component of Z, and C £ intersects 
some V C {U,W). 

Lemma 5.4. Z = lim £1)£2 Z. 

Proof. Denote e = (E2,£i). As lim e is a ring homomorphism from K(e)& 
to K, certainly lim £ Z C. Z. We shall show the reverse inclusion. Let x £ Z. 
We find a point x £ Z satisfying lim e a; = x. Note that p(Q(x)) G and 
p(Q(rr)) — £1 < 0, as £1 > e 2 > 0. On the other hand, as p(Q(X)) is not 
identically 0, for any < r G K the ball of radius r around x in K n contains 
a point y such that p{Q(y)) > 0. As x lies in the closure of the semialgebraic 
set F + defined by p(Q(X)) > 0, there exists a semialgebraic path 7 : [0, 1] — > 
K n such that 7(0) = x and 7((0, 1]) C F + , cf. Curve selection lemma 
Thm. 3.19]. As the image of a closed and bounded semialgebraic set under a 
continuous semialgebraic function on it is bounded, cf. Jj Thm. 3.20], 7QO, 1]) 
is bounded over K. 

Let 7 denote the extension of 7 to K(e). Then by jSJ Prop. 2.84] the set 
7Q0, 1]) is bounded over K(e). By the semialgebraic intermediate value theorem 
Prop. 3.4] the set X(r ) of all r G (0, r ) C K(e) satisfying p(Q(^(r))) = is a 
nonempty closed semialgebraic set. Choose r in the closest to interval of T(tq). 
Then lim £ r = 0, as To is arbitrary close to 0. As lim £ is a ring homomorphism, 
we have p((J(linif t"(t))) = 0. 

It remains to show that x = lim e 7(r). Identify r with the correspond- 
ing (representative of the) germ of semialgebraic continuous functions on K >0 
and think of j(t) as of composition 7 o r. To complete the proof, apply 
Lemma 3.21] that states that in this setting (7 a semialgebraic continuous 
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function on a closed bounded semialgebraic set over K and r an element of the 
extension of this set to one has 7(lim e r) = lim e (7 or). □ 

At this point we obtained a set of univariate representations u(eq,T) 6 
D[e ,T] n+3 (see Ql-ip ) for points in each connected component of Z. Now we 
get rid of the infinitesimal £q. Remove from u the polynomial gi responsible 
for X -coordinate, that is no longer needed, and apply [5, Alg. 12.46] (Removal 
of Infinitesimals), that consists of two steps. 

The first step is running the restricted elimination algorithm Alg. 12.43] 
with the input consisting of the polynomial / and the following polynomials: 

(5.5) f, f", . . . , /(**(/)-!), So CS(p(Q)) P(Q(-, • • • , —))■ 

go go 

It outputs a finite set S C D[e ] such that the degree of /, the number of roots 
of / in K, the number of common roots of / and h in IK, and the signs of h at 
the roots of /, for all h in (|5.5|) . are fixed on each connected component of the 
realization of any (realizable) sign condition on S. 
The second step computes, for each polynomial 

h(£ ) = h e £ e + h e+1 £ e +1 + ■■■ + h e+OJ £ e +u} e5, h t ^ 0, 

the Cauchy lower bound ^ ^4 | on the absolute value of its nonzero roots (see 
|5J Lemma 10.3]) and substitutes the minimum |, for a,b G D, of these bounds 
for Eq. The following remain unchanged upon substituting £ = | : 

o the number of real roots of /(|,T) and their Thorn encodings; 

o the signs of the polynomials in (|5.5J) at these roots. In particular the 
Thorn encodings of these roots will remain the same. 

By construction, the set of points represented by the it(|, T)'s intersects each 
connected component of Z. 

The number of arithmetic operations in D for these two steps is (dn)°( k \ 
according to jSI p. 462]. 

It remains to convert the T) = (/, g , g 2 , . . . , g n+ i)'s into real univariate 
representations by computing the Thorn encodings a for each real root of / 
using |5, Alg. 10.64] and Rem. 10.66]. The complexity of this procedure 
is (dn)°( k \ and for the case D = 7L the bitsizes of the intermediate data are 
bounded by (dn)°( k > times the bitsize of the input, by [loc.cit.]. 

This completes the proof of Theorem 11.21 



32 Grigoriev & Pasechnik 



Acknowledgements 

The second author is supported by NWO Grant 613.000.214. Parts of the work 
were completed while he held a position at THI/FB20 Informatik, University of 
Frankfurt (Germany) supported by DFG Grant SCHN-503/2-1. As well, he was 
partially supported by Universite de Rennes I (France) and the Mathematical 
Sciences Research Institute (USA). 

The authors thank Saugata Basu, Marie- Frangoise Roy, Nikolaj Vorobjov, 
and two anonymous referees for useful remarks, and Michael Kettner for care- 
fully reading a preliminary version of the text and tireless help in improving it. 
Sergey Fomin kindly provided references in Remark 12. 151 

References 

[1] M.-E. Alonso, E. Becker, M.-F. Roy, and T. Wormann. Zeros, multi- 
plicities, and idempotents for zero-dimensional systems. In Algorithms 
in algebraic geometry and applications (Santander, 1994), volume 143 of 
Progr. Math., pages 1-15. Birkhauser, Basel, 1996. 

[2] A. I. Barvinok. Feasibility testing for systems of real quadratic equations. 
Discrete Comput. Geom., 10(1):1— 13, 1993. 

[3] A. I. Barvinok. On the Betti numbers of semialgebraic sets defined by few 
quadratic inequalities. Math. Z., 225(2):231-244, 1997. 

[4] S. Basu, R. Pollack, and M.-F. Roy. On the combinatorial and algebraic 
complexity of quantifier elimination. J. ACM, 43(6)1002-1045, 1996. 

[5] S. Basu, R. Pollack, and M.-F. Roy. Algorithms in Real Algebraic Geom- 
etry. Springer- Verlag, 2003. 

[6] J. Bochnak, M. Coste, and M.-F. Roy. Real algebraic geometry. Springer- 
Verlag, Berlin, 1998. Translated from the 1987 French original, Revised 
by the authors. 

[7] D. Eisenbud. Commutative Algebra with a View Toward Algebraic Geometry. 
Springer- Verlag, New York, 1995. 

[8] D. Grigoriev. Complexity of deciding Tarski algebra. J. Symbolic Comput., 
5(l-2):65-108, 1988. 



Computing over quadratic maps I 33 



[9] D. Grigoriev and N. Vorobjov. Complexity lower bounds for computation 
trees with elementary transcendental function gates. Theoret. Comput. 
Sci, 157(2):185-214, 1996. 

[10] D. Grigoriev and N. N. Vorobjov, Jr. Solving systems of polynomial in- 
equalities in subexponential time. J. Symbolic Comput., 5 (1-2): 37-64, 
1988. 

[11] L. E. Heindel. Integer arithmetic algorithms for polynomial real zero de- 
termination. J. Assoc. Comput. Mach., 18:533-548, 1971. 

[12] U. J. J. Le Verrier. Sur les variations seculaires des elements elliptiques des 
sept planetes principales: Mercure, venus, la terre, mars, jupiter, saturne 
et uranus. J. Math. Pures Appl, 5:220-254, 1840. 

[13] I. G. Macdonald. Symmetric functions and Hall polynomials. Ox- 
ford Mathematical Monographs. The Clarendon Press, Oxford University 
Press, New York, second edition, 1995. With contributions by A. Zelevin- 
sky, Oxford Science Publications. 

[14] E. W. Mayr and A. R. Meyer. The complexity of the word problems for 
commutative semigroups and polynomial ideals. Adv. in Math., 46(3):305- 
329, 1982. 

[15] J. Renegar. On the computational complexity and geometry of the first- 
order theory of the reals. I- III. J. Symbolic Comput., 13(3):255-352, 1992. 

[16] F. Rouillier, M.-F. Roy, and M. Safey El Din. Finding at least one point 
in each connected component of a real algebraic set defined by a single 
equation. J. Complexity, 16:716-750, 2000. 

[17] J. von zur Gathen and J. Gerhard. Modern computer algebra. Cambridge 
University Press, New York, 1999. 



Manuscript received March 2, 2004 



34 Grigoriev & Pasechnik 



Dima Grigoriev 
IRMAR 

Universite de Rennes I 
Campus de Beaulieu 
35042 Rennes cedex 
France 

dimaOmath.univ-rennesl .fr 

http : //name .math. univ-rennesl . f r/dimitri . 
grigoriev/ 



Dmitrii V. Pasechnik 

Dept. E & OR and CentER 

Tilburg University 

P.O. Box 90153 

5000 LE Tilburg 

The Netherlands 

d . v . pasechnik@uvt . nl 

http : //center .uvt .nl/staff /pasechnik/ 



